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In this work, we investigate an exactly solvable two-leg spin ladder with three-spin interactions. 
We obtain analytically the finite-size corrections of the low-lying energies and determine the central 
charge as well as the scaling dimensions. The model considered in this work is in the same universality 
class of critical behavior of the XX chain with central charge c = 1. By using the correlation matrix 
method, we also study the finite-size corrections of the Renyi entropy of the ground state and of the 
excited states. Our results are in agreement with the predictions of the conformal field theory. 
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I. INTRODUCTION 

Perhaps one of the most popular and precise ways 
to determine the critical behavior of one-dimensional 
quantum systems is through the finite-size corrections 
of the energy. The machinery of conformal field the- 
ory (CFT) establishes that the ground state energy of 
a one-dimensional critical system of size L, under peri- 
odic boundary condition (PBC), behaves asymptotically 
asi^ 



Eo 
L 



6£ 2 



+ o(L- 2 ), 



(1) 



where v s is the sound velocity, is the bulk ground 
state energy per site, and c is the central charge. 

The mass gap amplitudes of the finite-size corrections 
of the higher energy states, are related with the scaling 
dimensions dP . There are, for each primary operator Op 
(J3 = 1, 2, ...) in the CFT, a tower of states E?., (L) in the 
spectrum of the Hamiltonian with asymptotic behavior^ 
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2ttv s 



{d p +j+f) + o(L- 1 ), (2) 



where = 0,1,2... . The above relations were used 
systematically to determine the universality class of crit- 
ical behavior of several models with great success. The 
success of these relations resides in the fact that it is 
possible to obtain the critical exponents, which are as- 
sociated with physical quantities in the thermodynamic 
limit, from finite systems. 

Recently, a great deal of excitement has been gen- 
erated due to the fact that the universality classes of 
critical behavior of one-dimensional systems can also be 
inferred from the finite-size corrections of the entangle- 
ment entropy of the ground state as well as of the excited 
states^— Below, we briefly report some important results 
about these corrections that we will explore in this work. 

Consider a quantum chain with L sites, described by 
a pure state whose density operator is p. Let us consider 



that the system is composed by the subsystems A with £ 
sites (£ = 1, . . . , L — 1) and B with L — £ sites. The Renyi 
entropy is defined as 



S a (L,£) 



1 



In Tr(p° 



(3) 



where pa — Trgp is the reduced density matrix of the 
subsystem A. The von Neumann entropy, also known 
as entanglement entropy, is the particular case a = 1. 
In the scaling regime 1 << £ << L, it is expected that 
for the critical one-dimensional systems with PBC, the 
Renyi entropy of the ground state behaves as 



S a (L,£) = S^ 1 (L,£) + S° a sc (L,£). 



(4) 



The first term in this equation, is the CFT prediction 
and is given by&l&~— 
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where c a is a non-universal constant. The second term 
is given by^r— 
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[aiding + g a {l - 5i,a) cos(kI + <f>)] 
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where p a is a new critical exponent, and a\ and g a are 
non-universal constants. The wave vector n and the 
phase <j) depend on the model. For instance, for the Ising 
model k = = 0, and for the spin-s XXZ chains at 
zero magnetic field n = it and = 0. For the systems 
with PBC, it is not expected oscillations in the von Neu- 
mann entropy iSi i 13 ' 16 The origin of the exponent p a in 
the above equation are the conical space-time singulari- 
ties produced in the conformal mapping used to describe 
the reduced density matrix p^ = Tr^p in the CFT ^2. 
This exponent is related to the scaling dimension x con of 



an operator (unknown) of the underlying CFT as follows 
p a = 2x con I aJ^- There are evidences that the exponent 
x con _ x ^ £ Qr q, > x, where x e is the scaling dimension of 
the energy operator, and that p\ = v = 2 for all models*^ 
On the other hand, the origin of the oscillating factor 
[cos(k£ + (f>)] is not yet completely understood, however 
it has been observed in systems with anti-ferromagnetic 
tendencies. 

The CFT prediction for the finite-size corrections of 
the Renyi entropy above were verified by several au- 
thors. Both the central charge ^ 15 ' 18 " — and the expo- 
nent p d 13 i 15 i 16 i 22 i 23 i 26 were obtained with amazing preci- 
sion. 

The Renyi entropies of the excited states also present 
very interesting universal behaviori 9 ' 24 ' 25 The a— Renyi 
entropy Sj, associated with an excitation of a primary 
operator T with conformal weight h, is related with the 
2a— point correlator of the operator T and its conjugate 

bj&24 



where 5£ s is the a— Renyi entropy of the ground state 
and 
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Figure 1: (Color online), (a) Schematic representation of the 
two- leg spin ladder. Ji, J2 and J3 are the couplings between 
the spins along the rung, the diagonal and the plaquettes, re- 
spectively. We considered a bipartite system as reported in 
the figure. The system of L rungs is divided into two subsys- 
tems of sizes £ and (L — £). (b) The spinless Hamiltonian 
The lines represent the hopping terms. 



4 a) = a- 2a ( h+ V 



U;ZoT[2nj/a}T^27T(j+£/L)/a} 
(T[0]Tt[27r£/L]) Q 



In particular, = 1 for the operators associated with 
compact excitations,^ 2 ^ and the Renyi entropy of these 
excitations are the same of the ground state, i. e. , 
5 T = S a s . On the other hand, for non-compact excita- 
tions the universal function depends on the partic- 
ular primary operator T. For instance, for the primary 
operator T = idej) of a free boson theory the universal 
function F^i was calculated analytically and is given 



9,24 
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where s — sin(7r£/2L). 

In this work, we intend to explore the above relations 
in order to investigate a two-leg spin ladder model, which 
is exactly soluble. The paper is organized as follows. In 
Sec. II, we present the model and its ground state phase 
diagram. In Sec. Ill, we determine the Renyi entropy 
of the model by using the correlation matrix method. 
Finally, in Sec. IV we summarize our results. 



II. PHASE DIAGRAM AND CRITICAL 
PROPERTIES OF A TWO-LEG SPIN MODEL 

We consider the following Hamiltonian defined on a 
two-leg geometry 
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where Jj(j = 1,2,3) are coupling constants [see Fig. 



QJa)], s \ = x iUi z ) are the spin-i operators at leg 
A = 1,2 and rung j, h is the magnetic field, and L 
is considered even. We investigate the above model 



with PBC, i. e., s 





X,j+L 



3 P 



and 



5? 



Note that this model has three-spin interactions. Models 



with multi-spin interactions like YYiLj+i s i s m+i (with 
P, j3' = x, y) can be mapped into a fcrmionic quadratic 
form, as first noted by Suzukij^l In the particular case 
of three-spin interactions, several variants of the Hamil- 
tonianjfS]) were considered in the literature^— (see also 
Ref. I35T ) . It is quite interesting to mention that spin-i 
Hamiltonians with three-spin interactions can be gener- 
ated using optical lattices^ 

The authors of Ref. HH also investigated the same 
model considered in this work. However, these authors 
focused in the effect of the lattice distortion in the in- 
finity system. Here, our emphasis is in the finite-size 
corrections (i) of the entanglement entropy and (ii) of 
the low-lying energy states in the critical regions. As we 
already mentioned, these corrections are related with the 
central charge and the scaling dimensions. 

For the sake of clarity, we briefly describe the 
procedure to diagonalize the Hamiltonian ([5]) exactly. 
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Figure 2: (Color online). Figures (a), (b) and (c) display the 
three typical profiles of the band dispersion, for h = 0, that 
appear in the regions I, II and III, respectively [see Fig. @]. 
kp are the Fermi momenta in the branches a = ± (see text). 



First, we use a Jordan- Wigncr transformation [c m = 
(s^ — isl^) Ylj <m {~2sj)} in order to map the Hamilto- 
nian in a spinless fermion chain with nearest-neighbor 
and next-neighbor hoppings [see Fig. Hfb)]. By using this 
transformation we obtain 
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H = J2 



Jicl. 



2L-2 T 

2 

3=1 



Sc. 



we can express the Hamiltonian (0) in the following di- 
agonal form 



H 



cr=± {fe} 



r (k)dpdl + hL, 



(11) 



where d^ 



j= 2 \a k ± 



Ji + J2 exp( — ik) 



b k , k = k(j) = 



•y/j 1 2 + ./|+2JiJ 2 COsfc 

with = (tt) if 7V> is odd (even) and j = 
0, ±1, ±(L - 2)/2, -i/2. The band dispersion is sep- 
arated in two branches (a = ±) given by 
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J? 



4 



2JiJ 2 cosfc. (12) 



Let us focus in the case of a zero magnetic field. After 
an exhaustive investigation of the band dispersion above 
with h — 0, we note that depending on the coupling con- 
stants Ji, J 2: and J3, the profiles of the band dispersion 
present three distinct behaviors. The main characteris- 
tics of these profiles are depicted in Figs. [2Ja)-(c). We 
label the coupling regions by regions I, II and III accord- 
ing to the Fermi momenta in the thermodynamic limit. 
In the region I (II) we have four (two) Fermi momenta, 
while in the region III the Fermi momenta always appear 
at kp = ±(ir — tt/L). We can observe from Fig. [2] that 
the regions I and II are gapless, whereas the region III is 
gapped. In the latter case, the lowest excitation energy 
appears in the sector s\ otal = — 1. The spin gap in this 
region is given by A = | Ji — J 2 \ — J3 ^ Fig. [3] displays 
these regions in parameter space J2/J3 vs J\/ 

It is interesting also to note that depending on the 
coupling constants the magnetization of the ground state 
may not be zero in regions I and II, differently from the 
region III where the magnetization is always zero if h = 0. 
The magnetization per site is given by jj- = ^ 5^ • < 

k p+ k F — 1/2 where the Fermi momenta 



gs\J2 s j\9s >- 

k F are determined solving the equations e ± (k 
the thermodynamic limit, and are given by 



±/u±> 



0. 



J1J2 - hJ 3 



± 



hJ\f F + hL, 



where Afp ~ ^2i=i c \ c i IS ^ ne particle number operator. 
Note that the original two-leg ladder system of size L is 
mapped into a chain with 2L sites. 

Finally, by using the following Fourier transformations 
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From the above equations we can see that for h = the 
(9) magnetization is not zero, in general. In Fig. [3J we also 
report the values of the intensity of magnetization m in 
the parameter space. 

Let us now investigate the critical regions I and II in 
more detail. The ground state energy of the Hamiltonian 
([TT]) is obtained adding particles in the energy levels be- 
low zero, i. e., 

(10) E -hL= J2 £ ~( fc )+ E e+ ( fc )> ( 14 ) 
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Figure 3: (Color online). The ground state phase diagram 
for h = 0. The regions I, II are gapless, whereas the region 
III is gapped. The right panel indicates the intensity of the 
magnetization per site m. 



where the last sum in this equation contributes only for 
the region I. We consider Nf = N t 



Np even and 
kp = (Afp — 1)tt/L are the Fermi momenta, where Afp 



is the particle number in the branch a . 

In order to determine the central charge c, we need 
to evaluate the finite-size corrections of the ground state 
energy [Eq. (fl"4)) ]. For a fixed value of the density p = ^ 
the leading finite-size correction can be obtained by using 
the Euler-Maclaurin formula^ which gives 
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If we compare Eqs. ([TJ and (fTS"]). we identify the central 
charge as c — 1 for each gapless mode. This value of c 
implies that the universality class of critical behavior of 
the Hamiltonian ([8]) is the same of the well known XX 
chain. Note that the critical behavior is described by 
two non-interacting conformal theories associated with 
the two branches a = ±. 

The excited states are obtained adding particles (holes) 
above (below) the Fermi level. Let us consider the fol- 
lowing excitations: (i) we can add (remove) I = 1+ + 1- 
particles above (below) the Fermi level, where I a [a = ±) 



is the number of particles added (removed) in the branch 
<r; and/or (ii) we can remove Q a particles in highest oc- 
cupied level of the left/right Fermi point of the branch 
a and add these particles in the lowest free levels of the 
right/left Fermi point of the same branch (see Table I). 
The finite-size corrections of these states can also be ob- 
tained by using the Euler-Maclaurin formula, as done for 
the ground state energy. The result obtained is the fol- 
lowing 
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where 5 a = mod [ mod (I a , 2) + mod (I, 2), 2] in the 
region I and 5 a = in the region II (note that in this 
region Vp = 0). We observe, again, from the above equa- 
tion that the Hamiltonian ([8]) is described by two non- 
interacting conformal theories. Similar expressions also 
have been found in models with more than one gapless 
mode^2r— Comparing Eqs. and ([XT]) we identify the 
scaling dimensions as 
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(18) 



The scaling dimensions above have a similar structure of 
the one of the Gaussian models Note that g?q x = x e = 1 
corresponds to the scaling dimension of the energy oper- 
ator, as in the XX chain. 



III. ENTANGLEMENT ENTROPY 

We are interested, now, in the finite-size corrections of 
the Rcnyi entropy, in the critical regions I and II. The 
leading finite-size corrections of Renyi entropy is related 
with the central charge, while the sub-leading corrections 
are governed by the exponent p a Below, we determine 
the central charge c and the exponent p a by using the 
universal behavior of these corrections. 

In the case of quadratic Hamiltonians, as the one con- 
sidered in this work, the Renyi entropy can be eval- 
uated through the correlation matrix metho d 44 ' 45 (see 
also Refs. and l47l). The main idea of this method 
is that the eigenvalues {vj} of the correlation matrix 
C p .q =< cjpCq >, p,q — 1, £; are related with the eigen- 
values of the reduced density matrix p^ . 44 ' 45 Due to this 
fact, the Renyi entropy can be expressed in terms of the 
eigenvalues {vj} in the following form 



S a (L,£) 
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1 



a 



5> [< + (i-^n 



(19) 



We will present data only for £ even due to the geometry 
of the two-leg ladder, although the odd sites also present 
the universal behavior predicted by the CFT. 
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The elements of the correlation matrix C Pi9 associated 
with the ground state wave function, can be easily eval- 
uated, and are given by 
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Similar expressions can also be obtained for the excited 
states. 

In the next two subsections, we determine numerically 
the eigenvalues of the correlation matrix C Viq in order to 
obtain the Rcnyi entropy of the ground state as well as 
of the excited states. 



A. Ground State 

Let us consider first the Renyi entropy of the ground 
state. In Fig. Sta), we present the von Neumann entropy 
Si, as function of I for systems of size L = 800, h = 
0, and three sets of couplings (Ji, J2,^)- These sets 
of couplings illustrate the behavior of S± in the three 
distinct regions of the phase diagram depicted in Fig. [3] 
As we observed in this figure, in the gapped phase (region 
III) the von Neumann entropy tends to a constant, as 
expected. On the other hand, for the critical regions I and 
II, Si increases in agreement with the CFT prediction 
[Eq. ©]. The black/blue lines in Fig. gja) are fits to 
our data using Eq. ([5]). The central charge c obtained 
through this fit is c = 1.00001 (c eff = 2.00021) for the 
sets of couplings associated with the critical region II 
(I). Similar results were acquired for several other sets 
of coupling for these two distinct critical regions. These 
results show that for the region II, the low-energy physics 
is describe by a CFT with central charge c = 1 , as we have 
already predicted [see Eq. (fT5|) ] . It is very interesting 
to note that in the region I, where we have two gapless 
modes (instead of one as in the region II), we obtain an 
effective central charge c e // = 2. This result resembles to 
the one found for the finite-size corrections of the ground 
state energy [Eq. (| 15[) ] . However, unlike the finite-size 
corrections of the energies, the entropy does not depend 
on the sound velocities. Then, it is expected that each 
gapless mode contributes for the finite-size corrections of 
the entropy as Eq. (|5|). For this reason, we get c e // = 
2c = 2. 

Another interesting feature of the von Neumann en- 
tropy in the region I, is that it presents an unexpected 
oscillation. In order to observe better these oscillations 
is convenient to define the following difference 
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Figure 4: (Color online), (a) The von Neumann entropy 
Si (L, I) of the ground state for systems of size L = 800, 
h = 0, Ji = 1 and some values of Ji and J 2 (see legend). 
The couplings associate with the black, blue, and red data 
correspond to the couplings belonging to the regions I, II and 
III, respectively [see Fig. ©]. (b) Results of the difference 
Di(L,£) for the same parameters of figure (a). Only few sites 
are presented. The solid lines in these figures connect the 
fitted points (see text). Insets: Show Si [figure (a)] and Di 
[figure (b)] for all sites. In order to show all data in the fig- 
ures, we added some constants in the values of Si and D\ , as 
indicated by the arrows. 
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where c e // =1(2) in the region II(I). From Eqs. (|3j)- (J6j) , 
it is expected that the difference D a behaves as 
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In Fig. QJb), we show D\ (L, t) for the two sets of coupling 
parameters presented in Fig. |3Ja). If Ji ^ J 2 , oscilla- 
tions between the odd and even sites appear naturally in 
the entropy. The oscillation observed in Fig. |U[b), for the 
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set of coupling belonging to the Region I, is not related 
with the "dimerization", is present even for J x = J 2 . Ac- 
cording to a previous conjecture of Xavier and Alcaraz,— 
based in studies of several models with just one gapless 
mode, the von Neumann entropy of single interval should 
not present these oscillations. Our present results, show 
that the same conjecture does not apply for models with 
more than one gapless mode. These authors also conjec- 
tured that the sub-leading finite-size correction of the von 
Neumann entropy decays with the exponent p\ = v = 2. 
Indeed, we also have observed this decay in the region II. 
The blue solid line in Fig. |U[b) is the fit to our data as- 
suming that Di(L, £) behaves as Eq. (|2Tj) . The exponent 
that we get by this fit is p\ = 2.0001, in agreement with 
the previous conjecture. 

In the region I, where we have two gapless mode, we 
fit our data assuming that Di(L,£) has an oscillating 
term different of Eq. ([5]). In this region, we had to re- 
place gi(l — $i, a ) cos (n£ + <f>) in Eq. © by the following 

term: gi cos (£kp) + g^ cos (£kp^j + g[ 3 ^ cos ( — ~^ kp ' 
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in order to fit our data. By using this 

oscillating term, we are able to fit perfectly our data, 
as shown in Fig. EJb). The exponent obtained in this 
later case is p\ = v = 2.0007. Similar results were also 
observed for other couplings. The motivation for us con- 
sider the above oscillating term is given below. These re- 
sults indicate that the sub-leading correction of the von 
Neumann entropy indeed has the exponent p\ = v = 2 
even for models with more than one gapless modes. 

Let us now consider the a— Renyi entropies of the 
ground state with a > 1. As illustration, we present 
in Fig. [5£a), S a (L,£) as function of £ for systems of size 
L = 800 and L = 40, two values of a, and some cou- 
plings. In the case that a > 1, it is expected oscillations 
in the entropy even for systems with PBC^ as men- 
tioned in the introduction. Moreover, it is expected that 
the amplitude of the oscillations decrease, as the system 
size increases [see Eq. ([6])]. Indeed, this feature is clearly 
noted in a system with smaller size, as depicted in the 
inset of Fig. [5ja) . 

Once again, in order to observe better these oscillations 
we investigate the difference D a (L,£). In Fig. 03h), we 
show D a (L,£) vs £ for the same data presented in Fig. 
[5^a) with a = 3. Let us first discuss the data of the Renyi 
entropy in the region II. In this region, we are able to fit 
our data by using Eq. © with k = kp and = 0. The 
exponent we get from this fit is pz = 0.666. Similar fits 
of S a (not shown) for a = 2 and a — 4 give p2 — 0.999 
and p4 = 0.497, respectively. These results indicate that 
p a = 2/ a, as in the XX chain. 

On the other hand, it is not possible to fit the data 
of the Renyi entropy, in the region I, if we consider that 
the oscillating term behaves just as cos(«<? + </>), as it was 
done in the region II. In order to get some insight about 
these oscillations in the region I, we calculate analytically 
the spin-spin correlation function (s2n s 2n+2i)i which is 
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Figure 5: (Color online), (a) The Renyi entropy of the ground 
state for systems of size L = 800, h = 0, J3 = 1, a = 3, and 
two values of Ji and J2 (see legend) . The couplings associate 
with the black and blue lines, correspond to couplings belong- 
ing to the regions I, and II, respectively [see Fig. (|3J]. Inset: 
Results for L = 40. (b) The difference D a (L,l) for the same 
parameters of figure (a). Only few sites are presented. The 
symbols in the figures (a) and (b) are the numerical data and 
the solid lines connect the fitted points (see text). 



given by 



S 2n s 2(n+t) 



AL 2 sin 2 {-K£/2L) 



kF+kF £)+co S 



1- - cos (£k F ) 



Motivated by the fact that it is expected that the os- 
cillations of the entropy are associated with the anti- 
ferromagnetic nature of the Hamiltonian^ we assume 
that the difference D a (L,£), in the region I, behaves as 
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Table I: Schematic representations of some states. The open 
(closed) cycles are the empty (occupied) levels. The two rows 
of cycles inside the symbol {} denotes the two branches, while 
the vertical lines show the positions of the Fermi momenta. 
We also present the excitations in terms of 7 CT and Q a [see Eq. 
I|17p] , and the notation used by the authors of Ref. [24]. In the 
later case {h\hi... : piP2...) a represents an excited state with 
holes (particles) in the his (pis) allowed momentum values 
below (above) the Fermi point of the branch a. 



Indeed, as illustrated in Fig. OJb) , we obtain a very nice 
fit of our data if we assume that D a behaves as the above 
equation. The exponent we get from the fit is p% = 0.667. 

Similar results, as the ones presented above, were ob- 
tained for several other values of m and coupling param- 
eters. In addition to that, the results for some values of 
a strongly support that the exponent p a = 2/ a, as hap- 
pens in the XX chainJ^ These results are in agreement 
with the conjecture^ that the exponent p a is related with 
the dimension of the energy operator x e by p a — 2x e /a, 
which is x e = 1 in the present model. 



B. Excited states 

Finally, we consider the Renyi entropy of the excited 
states. Differently of the Renyi entropy of the ground 
state, where several studies indeed confirmed the univer- 
sal behavior predicted by the CFT, few works considered 
the excited states ,2i24 ~ 26 i 48 ~ £li Certainly, more studies are 
highly desired in this front and in this subsection we in- 
tent to provided a study of the Renyi entropy of several 
excited states of the Hamiltonian ([5J . 

First, we consider the Renyi entropy associated with 
compact excitations. The compact excitation are those 
that do not present holes in the momentum space. In 
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Figure 6: (Color online), (a) AS2 vs £ for some particle-hole 
excitations for systems of size L = 1000, h = 0, and some 
couplings (see legend). The symbols are the numerical data 
and the solid lines connect the fitted data (see text). The 
solid curves from the top to down have values of n: 4,3,2 and 
1 [see Eq. [23] ■ The blue(red) symbols are data from the region 
I(II). Only few sites are presented. In order to show all data 
in the figure, we added some constants in the values of AS2, 
as indicated by the arrows. We also present the excitations in 
terms of the notation of Ref. [24|. Inset: AS2 for all sites, (b) 
A Si vs 2h_(^/2Z/) 2 for several excitations for systems of size 
L = 10000 (see Table II). 



Table I, we depict few examples of such excitations. We 
found that AS a {L,£) = S e a xct - Sg s < 1CT 3 for these 
compact excitations, as well as many others not presented 
in this table. These results are in agreement with the 
prediction, done by the authors of the Refs. HI and H3, 
that S^ xc = Sg s for the compact excitations. 

The Renyi entropies associated with the non-compact 
excitations are, in general, different from the ground 
state. Let us consider particlc-holc excitations in the 
right/left Fermi points of the branch a. In Fig. EJa), we 
present AS2 as function of £ for some these excitations 
for systems with 2L sites. As observed in this figure, 
the Renyi entropies associated with these non-compact 
excitations increase with I. Moreover our results, show 
that the AS2 data can be perfectly fitted (except by the 
unusual oscillations) if wc assume that [see Eq.(J7J)] 
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AS 2 = -nF^l = -nln (1 - 2s 2 



3s 4 - 2s 6 



8s a 



REGION II 
I- Q- i- 



(23) 

where s = sin(7r£/4L) and n is the number of particle- 
hole excitations associated with the excited state. 

In the case of just one gapless mode, the Eq. (|23|) was 
determined analytically for the free boson theory^ Note 
that for a free boson theory with two decoupled field </> CT 
(a = ±) the Eq. is valid - 

Finally, we consider ASi (L, £) = Sf xc — Sf s for several 
excitations. In particular we focus on the behavior of 
A 5i in the regime that (f/L) < 1. In this regime it is 
expected that ASi, f° r a system with 2L sites, behaves 



o[(£/L)% 



(24) 



where d@ is the scaling dimension of the operator asso- 
cociate with the excitation. This behavior of ASi above 
was proofed^ in the case of the primay operator id(j> and 
verified numerically for the descendant operators.— 

The scaling dimensions of the XX chain have the same 
structure of the Gaussian model 3 - (with compactification 
radius R = 2) and are given by di t Q t j t j> = xi.q + j + j', 

where xj.q = (-j- + Q 2 ) arc the scaling dimensions of 
the primary operators and j,j' = 0,1,2... [see Eq. fl2J)]. 
Note that the energy associated with an operator with 
dimension dj q j j> (in the Gaussian model) has momen- 
tum AP}£ = ^-{IQ+j - /)■- We also found a sim- 
ilar result. The momentum of a state associated with 
the scaling dimension J^ tT=± (dj q + j a + j' a ) is given 
by AP = Z<,±*PL QM , where AP£, Q ^,, = 
X [I<t(Q<7 + <W2) + ja — j' a ]- The momenta of the con- 
figurations presented in Table II are in perfect agreement 
with the values of AP. 

In Fig. fflb), we show ASi as function of ^-(£/2L) 2 
for several excitations, in the regime that (£/L) <C 1 for 
a system with 2L sites. The schematic representation 
of these non-compact excitations are also presented in 
terms of I a , Q a and j a in Table II. The solid lines in Fig. 
[6]Jb) correspond to Eq. (|24|) with the values of dP shown 
in parentheses. As observed in this figure, the values 
of S 3 that fit our data to Eq. are consistent with 

d® = J2a=±(j<? + j'a)- Slightly different than expected 
from the Eq. (|24[) . In order to better understand this 
discrepancy, we also reproduce the results of the same 
excited states considered in the Refs. and HH for the 
XX chain. Indeed, we also get the same values of d 13 
obtained in these references. For example, for the excita- 
tion associated with the operator idcf>, i. e. (1 : 1), we got 
d@ = 1. However, we did not interpret this value as the 
scaling dimension of the primary operator id(f> (which is 
1), due to the momentum of the excitation (1 : 1). From 
the CFT point of view this primary operator has momen- 
tum 0. However, the excitation (1 : 1) has momentum 
AP = 2ir/L. These results show that a deeper analysis 
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1 2-10 2 
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Table II: Schematic representations of the excited states con- 
sidered in Fig. [6jb). We also present the excitations in terms 
of To- and Q a [see Eq. (|17[)] as well as the notation of Ref. l24l . 



of the Renyi entropy of excited states associated with the 
descendent operators should be done in order to under- 
stand the exact dependence of ASi, in the regime that 
(£/L) <C 1. Below, we comment on two points that may 
help infer the exact form of the factor dP that appears in 
(f24"l) . First, note that Eq. (|24| was obtained exactly for 
the primary operators. For instance, the operator id(j> 
is primary in terms of the Virasoro algebra. However, 
for a larger algebra as the Kac-Moody algebra (which is 
the correct one for the XXZ chain) the operator id<j> can 
be seen as a descendent operator. So, the result of Eq. 
(|24l) independent whether the operator id(f> is primary or 
not. Second, our results suggest that the operators asso- 
ciated with configurations as the ones depicted in Table 
II may be express by a product of vertex operators with 
scaling dimensions di t Q times a product of derivatives of 
the fields <fi and <fi with scaling dimension j + j' , and this 
would lead a factor d 13 = j + j' in Eq. (pM)) ^ Further in- 
vestigations are needed in order to clarify the exact form 
of ASi in the limit (£/L) — > for a general excited state. 
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IV. CONCLUSIONS 

In this work, we considered a toy model with three 
spin-interactions defined on a two-leg geometry which is 
exactly solvable by means of a Jordan- Wigner transfor- 
mation. We were able to obtain analytically the finite- 
size corrections of the low-lying energies, and for this rea- 
son we get the exact values of the central charge and the 
scaling dimensions. The toy model investigated [Eq. ©] 
is in the same universality class of critical behavior of 
the XX chain, with central charge c = 1. Moreover, 
the scaling dimensions of the toy model present a similar 
structure as the one of the XX chain [see Eq. (TT51) ]. 

We also acquired the central charge by analyzing the 
finite-size corrections of the Renyi entropy. In the region 
of the phase diagram where there is one gapless mode, 
we get c = 1. Whereas in the region where there are 
two gapless modes, we found that Renyi entropy behaves 
asymptotically as 
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— sin 


:?)] 




OLj 




7T ' 





where c e // = 2c = 2. For systems with n g i gapless 
modes it is expected that c e ff = n g ic, at least for non- 
interaction systems (see also Fagotti in Ref. [34| for a simi- 



lar discussion). It is interesting to mention that an exten- 
sion of the present model for AT- leg ladders with (N + 1) 
spin interactions can also be done and would lead to an 
effective central charge c e // = TV (at least for some re- 
gions of the parameter space) . This result shows that the 
von Neumann entropy of the TV- leg spin ladders of size L 
and L gapless modes should behave as Si ~ L ln(£) + aL, 
this simply argument shows that entropic area law may 
be violated in two-dimensional critical systems. Indeed, 
violations of the entropic area law was observed in sys- 
tems in dimension higher than on o 53 ^ 54 

We also study the Renyi entropy of the excited states. 
Our results of the Renyi entropy associated with the com- 
pacts as well as the non-compact excitations are in agree- 
ment with the recent prediction done by the authors of 
Refs. landH. 
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